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INTRODUCTION 

The object of this paper is to investigate the deviation of the solution of the 
differential equation 


y' s f( t, y), y(0) = y 0 


from the solution of the perturbed system 


z' = f(t, z) + g(t, z), z(0) = ZQ 


where f and g are continuous functions on I x R n into R n , where 
I = (o, a) or I = (o, °°). These functions are assumed to satisfy the Lipschitz 
condition in the variable z. 

The space Lip(I) of all such functions with suitable norms forms a Banach 
space. Introducing a suitable norm in the space of continuous functions C(I), 
one can reduce the problem to an equivalent problem in terminology of operators 
in such spaoes. 

For the sake of completeness of the presentation, a theorem on existence and 
uniqueness of the solution is presented by means of Banach space technique. 

Norm estimates on the rate of growth of such solutions are found. As a conse- 
quence, estimates of deviation of a solution due to perturbation are obtained. 
Continuity of the solution on the initial data and on the perturbation is 
established. 

As an illustration a nonlinear perturbation of the harmonic oscillator is 
considered. Also, a perturbation of equations of the restricted three-body 
problem linearized at libration point is considered. 


DISCUSSION 

Let Y be the R n space and | | a norm in it. Readers familiar with vectoral 
integration in Banach spaces (ref. 1) may assume that Y is a Banach space. 

Let I be a closed interval of the form (o,a), where a > 0, or (o, 00 ). Let 
f be a continuous function from the product I x Y into Y such that 

| f(t, y) - f(t, z)| < l| y - z| 

for all y, z e Y and all tel, and moreover, 
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I f(t, 0)| < L 


for all tel. 

Suoh functions will be oalled Lipschitzian in the sequel. Let 1 1 f 1 1 denote 
the infinum of all constants L appearing in the previous inequalities. De- 
note by Lip(I) the set of all such Lipschitzian functions from I x Y into Y, 
The pair Lip(l), |) J| represents a Banach space. 

The object of this paper will be to establish how much the solution of the equa- 
tion 


y'(t) = f(t, y(t) ) for tel 


satisfying the initial condition y(0) = yo will change in a suitable norm in 
the space of continuous functions if the function f is perturbed by a 
Lipschitzian function g of small norm. 

Let C k (I) denote the set of all continuous functions h from the interval 
I into the Banach space Y such that the number 

sup 1 1 h(t)| e" kfc : t e 1 1 = 1 1 h| | k 

is finite. Notice that the pair 0^(1), || | | ^ forms r Banach space for every 
k >_ 0. Moreover if the interval I is bounded; i.e., a< 00 , then all the 
spaces Cq(I) = 0^(1) are equal, and the norms || |)o and || ||^ are equiva- 
lent for any k. If the interval I is unbounded, then C^(I) is contained 
in C m (I) if k < m. 

For the sake of completeness of the presentation, the following theorem shall 
be proved. 

Theorem 1 . For every Lipschitzian function f e Lip(I) and every initial value 
yo £ Y there exists a constant kQ such that the differential equation 


y'(t) = f(t, y(t)) for all tel 


( 1 . 1 ) 


with the initial condition 


y(o) = yo 


( 1 . 2 ) 


has a unique solution in every space Cyr(I) such that k >_ kQ. 
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Proof . Using veotoral Integration, it can be established that the equation (1.1) 
with the initial condition (eq. 1.2) is equivalent to the integral equation 


t 

y(t) = yo + J f(s, y(s)) ds for all t e I (1.3) 

0 

Let F denote the operator defined by the right side of the equation and 
considered over the space C^d) ; l.e., h = F(y), where h and y are funo 
tions and y e C^d), if and onl y if » 


t 

h(t) = yo + f f(s, y(s))ds for all t e I 
0 


It can be proved that the function h is continuous for every function y e 0^(1). 
Lemma . The operator F maps the space C^d) into itself and 

| |F(y) - f(z) || k < |K|| | |y “ z | |k< 

for all y, z e 0^(1) and all k > 0. 

Proof of the Lemma . Take any y, z e 0^(1) and let y = F(y) and z = F(z). 

Then by definition of the operator F 

t 

|y(t) - z( t) | = | J (f(s, y(s)) - f ( s , z(s)))ds| 

0 


It follows from the definition of the space Lip (I) that 


|f(s, y(s)) - f(s, z(s))| ^ ||f|| |y(s) - z(s)| 


for all s e I* From the definition of the norm ^ we get 


3 


80 FM 21 


| y(s) - z(s)| <, 1 1 y - z| | k e ks for all 3 e I. 
Thua we get the estimate 

t 

| y(t) - z(t)| < 1 1 f| | 1 1 y - z| | k J e ks ds 

0 


for all t e X, Became 
t 

J e ks ds s (e kt - 1)/k <. e kt 7 k 
0 

for all t e I, we get 


e -kt| y(t) - z (t)| < 1 1 f 1 1 1 1 y - z| | k/k 


for all tel. Taking the supremum over all t. e X of the left aide of this 
inequality, one gets 

1 1 y - z| I k < 1 1 f| I 1 1 y - A I k/K 


and this proves the lemma. 

Now let kg = 1 1 f | | . Then for k > kg the operator F is a contraction of 
the Banach space Ck(I). Thus, from Banach's fixed point theorem (ref. 1 ), 
we get the existence of a unique point ye 0^(1) such that 


y = F(y) 


This proves the theorem. 
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Theorem 2 . Let y, a be two funotlons from the space 0^(1) and c > 0. If 


t 

|y(t)| < I a< t) | + o J |y(s)|ds for all tel (2,1) 

0 


then 


1 1 y I Ik 5 I Ul Ik + (o/k) 1 1 y I Ik for k > 0 


( 2 . 2 ) 


Proof . From the definition of the norm || ll^ it follows that 


t t 

/ |y(s)|ds< f 1 1 y 1 1 k eks ds < ||y|lke kfc /k 

0 0 

for all t e I and k > 0. 

Thus, from inequality (eq. 2.1) we have 

I y ( t) | < |a(t)| + o| | y { lke kfc /k for all tel 

Multiplying both sides by e" k ^ and using the definition of the norm || ||k 
we get 

e-kt|y( t )| < ( | a | | k + (c/k)||y||k 

for all b e I. 

Hence 


|y| Ik 1 I Ul Ik + (o/k)| | y | Ik 


Theorem 3 . Let f be a function from the space Lip(I), yg e Y, and let 
y e Ck(I), where k > ||f||, be the solution of the initial value problem 
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y'(t) s f(t, y(t) ) for all tel 
y(0) = yo 

Then 

l|y||k < <! yol + 1 1 f 1 1 /k ) / ( 1 - ||f| J/k) 

In other words, 

|y(t)| < e kfc ( | y 0 | + ||f||/k)/(1 - 1 1 f| | /k) 

for tel. 

Proof . 

The following estimate is derived from the triangle inequality and from 
the definition of the space Lip(I): 

| r(t, y) | < |f(t, y) - f(t, o)| + |f(t, o)| 

< l|f|| |y - o| + ||f|| 

< || f I| ( ( y| + 1) for all tel, y e y. 

Thus, for any y e C^d) satisfying the initial value problem, 

t 

y(t) = yo + f f(s, y(s))ds for all tel 
0 

This yields 

t 

|y(t)| < | yo | + j J f(s, y(s))dsj for all tel 

0 

From the previous estimate we get 
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f(a, y(a))| < I ! f| I (I y(a)| + 1) for ail a e I 


Thus 


t 

I y(t)| < (| yol + 1 1 r| I t) + 1 1 f| I J I y(a)| do 

0 


Using Theorem 2 and noticing that the norm of the function a, given by the for- 
mula 


a(t) = |y 0 | + || f J | t for all tel 


satisfies the estimate 


I M I k i I yoi + 1 1 f| I /k 


we get the required inequality. 

We shall prove the following theorem because of its simplicity. It is suffi- 
cient in the investigation of continuity of solution with respect to pertv rbatlon. 
Finer results; will be given in reference 2. 

Perturbation Theorem . Let f and g be two Lipschitzian functions from the 
Banach space Lip(I). Let yg, zq e y and let y, z e 0^(1), where 
k > ko - 1 1 f 1 1 + ) I g| | , be the solutions of the initial value problems 

y'(t) = f(t, y(t)) for all tel, y(0) = yp 

and 

z'(t) s f(t, z(t)) + g(t, z( t) ) for all tel, z(0) = zq 
Then the following inequality holds 

lly - zllk i Y 


where 
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Y = (I yo - -ol + 1 1 g! I o/k)/( i - (Ufll + llgll )/k) 
c s (lyol + 1)/(1 - 1 1 f| I /k) 
or equivalently 

I y(t) - z(t)| < Y e^ for all tU 

Proof . Using the integral equivalents of the above initial value problems we 
get 


t t 

y(U) s y 0 + j (f(s, y(s)) + g(s, y(s) )ds - J g(s, y(a))ds 
0 0 

t 

z(t) s zq + j (f(s, z(s) ) + g(s, z(s))ds 
0 

for all t e 1. 

Subtracting both sides of these equalities and using the triangle Inequality, 
one gets 


t 

| y(t) - z(t)| < | yo - zqI + J (I I f| I + 1 1 gl I )(l y(s) - z(s)| )ds 

0 

t 

+ / I Igl I (ly(s)l + 1)ds 
0 


for all t e X. 

The last term of the Inequality can be estimated as follows 
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t 

/ IUII <1 y(»)| + Ddfl i IUII 

0 


fc 

/ (| I y| I k° k8 + e k8 )^» 
0 


< lUIUIUIIk + 1)e kfc /k 


for all be X» 

Similarly one gets the estimate 


b 

/ <ll f ll + I j ej | ) |y(») - »<»)) d 8 < (||f|| + || S ||) 1 1 y - zjj k eUt/k 
0 


These estlmabes yield 

|y(b) - z(b)| < | yo - ®o| + (||f|| + | U| | )(||y - z|| k )e kt /k 
+ ||g||C||y||k + 1)e kfc /k 


for all b e X. 

Multiplying both sides of this inequality by e~ kfc and using the definition of 
the norm 1 1 1 1 k one gets 

||y - z||k < Uo - so| + <11*11 + 1 1 g| | )<| ) y - z| |k) /k 
+ I Ull (| |y| |k + U/k 

From theorem 3 we have the estimate 


1 1 y| Ik < <|yol + ||f||/ k )/d - lUH/to 


Using it and taking the term containing ||y - z from the right side onto 
the left side in the previous inequality, one gets,’ after simplifications, the 
following inequality 


|y - *||k < Y 
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where 

Y = (|yo - s=o I + I Is I |c/k)/( 1 - (||f|| + I Is 1 1 )/k) 

and 


c = ( |yo I + 1)/(1 - I |f | |/k) 


This proves the perturbation theorem. 

Corollary . Let g n e Lip(I) be a sequenoe such that ||g n || * 0, 

Nsnll 1 b for a11 n and 2 n e Y such that |z n - yo I * 0. Let h n e C^CI) 
be the solution of the perturbed initial value problem 

hn(t) = f(t, h n ( t) ) + g n (t, h n ( t) ) for all t e I 

h n (0) = z n 

Assume that y is the solution of y'(t) = f(t, y(t)) for tel with the 
initial condition y (0) = yg. Then 


I ly “ h n | Ik 0 


when n tends to infinity for any fixed k > ||f|| + b. 

This shows that the solution y e C^(I) of the initial value problem depends 
continuously on the initial value yg and Lipschitzian perturbations of the 
function f. 

Examples 

Nonlinear perturbation of harmonic oscillator 

Consider the following second-order differential equation 


x" + = e(cos x + sin x') 

x(0) = xg, x'(0) s Vq (A) 
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Treating the right side as a perturbation, find an estimate how much the solu- 
tion of the equation 


x" + wx = 0, x(0) = xq, x'(0) = v o 


(B) 


differs from the solution of the equation (A). 

Introducing the vector variable y = (yi, yz* b y the formulas yi = x, y 2 = 
and writing the equation in the normal form one gets the system 


y] = y2 

y^ = -w 2 yi + e(cos yi + sin >* 2 ) 
with initial conditions 


yi (0) = x 0 , y2(0) = vo 

Equivalently in vector form we get 


y' = f(y) +eg(y), y(o) = (x 0 , vo) 

where 

f(y) = (y 2 - a) 2 y 1 ) 

g(y) = (0, cos yi + sin yz) 
for all vectors y = (y^ y 2 ) e R 2 . 

The Lipschitz norm of the function f can be estimated as follows. Introduce 
in R 2 the norm 

I y| = I y2| + I y2| 

for all y = (yi , y 2 ) • 
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Let z = f(y) and z s f(y) then 

|z ~ z| = |z ■) ** zi | + |z 2 - z 2 | = I yg - y2 1 + w ^lyi - yi I 
- ^( w )(lyi - yi I + Iy2 - y2 I ) < ^( w ) |y - y| 

where X(W) = 0)2 if to > -j and \(w) - 1 If 0 < w < 1. 

This yields the estimate 

If(y) - f(y)l < X(W)|y - y\ 

for all y, ye r 2. Because f (0 ) = 0, we get for the Lipschitz norm of f 

| If | | < X(w) 

Similarly, let z = g(y) and z = g (y). 

Then 

|z — Z | = 1 55 -J — Z 1 I + |z2 “ Z2 I 

= 0 + loos y-| cos yi + sin y2 - sin y 2 I 

From the Chauchy theorem of calculus 

cos yi - cos = (-sin o-, ) (y n - yi ) 
sin y 2 - sin y2 = (cos c 2 )(y2 - yz) 

for some intermediate values ci and c 2 . 

These values give the estimates 

|cos yi - cos yi | < | yi - 7i I 
|sin y 2 - sin y 2 | < |y 2 - y 2 | 

Thus from (C) one gets 
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z - z] 1 1 yi - yi| + | y2 - ?2l s I y - yl 


i * 0 • f 

| g(y) - g(y)| 1 1 y - y| 

for all y, y e R 2 . 

Because 

|g(0)| = | (0, 1)| < 1 

we get the estimate for the Lipschitz norm of g 

Nell < 1 

Just from the perturbation theorem one gets the estimates 
Y <. (ec/k)/(1 - (1 + e)/k) 


where 

o < (| x 0 | + | vo| + 1 )/(1 - 1/k) 
if k > 1 + e . 

To be concrete let k = 2, e = 0.5, a) = 1 , Xq = 0, Vq = 1. Then 
Y > (0.5 • 4/2)/(1 - 1.5/2) = 4 

Thus 

|y(t) - z(t)| < 4e 2 ^ for all t > 0 
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The solution of the unperturbed system y' s f(y), y(0) s (0,1) Is given by 

y-j ( t) s sin (t) 
y 2 (t) = oos (t) 

for all t e k. Comparing it with the solution of the perturbed system 
z' = f(z) + gg(z), z( 0 ) = ( 0 , 1 ), one gets for the funotion 

h(t) = | y ( t) - z(t)| = (sin (t) - z-|(t)| + [cos (t) - zg(t)| 

the values in table I. 

Perturbation of restricted three-body problem linearized near libratlon point 
Consider the equations of restricted three-body problem 

x - 2y = Bx + Ay 
y - 2x = Ax + 33y 

and perturbed equations 

x - 2y = Bx + Ay + e k(x) 
y - 2x = Ax + 3By + e k(y) 

where k(u) = u/(7 + u 2 ) for all u e R. 

Introducing a new vector variable z = tz-| , Z 2 , Z 3 , Z 4 ) by the formulas 

x = zi, x = Z 2 , y = Z 3 , y = z^ and writing the equations in the normal 
form one gets 

z s f(z) + e g(z) 


where 

f(z) = z for all z e R 4 
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and 


4 


being the matrix 


I 



0 , 1 , 0 , 0 
B, 0, A, 2 
0 , 0 , 0 , 1 
A, -2, 3B, 0 


A 



2M) , 



y S 1/82.45 


and 


g(z) 


0 

k(z2) 

0 

k(z 3 ) 


for all z e R‘ 4 , Notice that f(0) = 0 and g(0) = 
mate of the Lipschitz norm 1 1 f 1 1 and 1 1 g| | it is 
mate of the Lipschitz constant for the functions f 

any z, z e R^ and estimate | f(2) - f(Z)| . 

Using the norm | z| = |zi| + | Z 2 I + | zg| + | zJ for 

mate ! 


| f ( z ) - f ( z ) | < (A + B) | z 1 - z 1 1 + 3| Z 2 - Z 2 I 
+ 3| z 4 - z 4 | 


Thus, if 


0. Thus, to find an esti- 
sufficient to find an esti 
and g. To this end take 


z e R 1 *, we get the esti 


+ (A + 3B) zg - zg| 


A = max {3, A + 3B} = A + 3B s 3*5175272 
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we get 

| f(z) - f(z)| < X | z - z| for all z, z e R 1 * 

This yields | | f | | < X . 

To estimate the Lipschitz norm of the function g let us first find the 
Lipschitz norm of the function k. Notice that 

Mx) = |k'(x)| = | (1 - x 2 ) | /( 1 + x 2 ) 2 

for all x. 

The function % has an absolute maximum on the interval (- 00 , 00 ) at x = 0 . 
Thus, 

| k» (x) | < MO) s 1 
for all x e (-00,00). 

Using the Cauchy intermediate value theorem, we get 
| k(y) - k(y)| = |k'(c)| |y - -I < |y - y| 
for all y, y e R. 

Now from the definition of the function g one gets 

| g( z) - g(z)| < |k(zg) - k(zg) | + | k( Z3) - k(z3)| 

< |z 2 - z 2 | + 1 23 - Z 3 | 

< I Z1 - Z1 | + | Z2 - Z2 1 + 1 23 - 23 1 + I z n - Zjj I 
= |z - z| 

for all z, z e R 1 *. This yields 


|g| | < 1 
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The estimates for the deviation I y ( t) - z(t)l one gets similarly as in the 
previous example. 

For formulas on Lipsohitz constant for the case of Euclidean norm and for finer 
error estimates, see reference 2. 


CONCLUDING REMARKS 

Though the estimates on deviation obtained in this paper are simple, the 
results they yield for longer periods of time will be very pessimistic. 
Finer estimates can be found in reference 2. That reference also contains 
general methods of computing Lipsohitz constants. 
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TABLE X.- COMPARISON OF COMPUTED GLOBAL ERROR WITH PREDICTED BOUND 


t 


h(t) 

! 

! 

! 

R exp(2fc) 

0.0 


0.0000000 

! 

! 

4.0000000 

.1 


.0976182 

! 

4.8856111 

.2 


.2050097 

1 

5.9672988 

• 3 


.3199916 

! 

7.2884752 



.4401596 

! 

8.9021636 

.5 


.5629649 

! 

10.8731271 

,6 


.6857933 

! 

13.2804674 

.7 


.8060412 

! 

16.2207992 

• 8 


.9211782 


19-8121288 

• 9 


1,0287934 

! 

24.1985884 

1.0 


1.1266198 

1 

29.5562222 

1.1 


1.2125354 

! 

36.1000504 

1.2 


1.2845469 

! 

44.0927005 

1.3 


1.3407616 

1 

53.8549452 

1.4 


1.3793622 

! 

65.7785778 

1.5 


1.3985967 

! 

80.3421354 

1.6 


1.3968026 

! 

98.1301050 

1.7 


1.3724767 

! 

119.8563795 

1.8 


1.3244038 

! 

146.3929100 

1.9 


1.2518431 

! 

178.8047028 

2.0 


1.1547580 

! 

218.3925552 

2.1 


1.0340512 

! 

266.7452583 

2.2 


1.0904912 

! 

325.8033905 

2.3 


1.2169999 

! 

397.9371529 

2.4 


1.3288932 

! 

486.0415306 

2.5 


1.4217193 

! 

593.6524582 

2.6 


1.4916894 

I 

725.0887451 

2.7 


1.5360930 

! 

885.6253815 

2.8 


1.5535694 

! 

1081.7052612 

2.9 


1.5441907 

! 

1321.1977844 

3.0 


1.5093642 

1 

1613.7145996 

3.1 


1.4515959 

! 

1970.9954376 

3.2 


1.3741666 

! 

2407.3792419 

3.3 


1.2807751 


2940.3796082 

3.4 


1.1751914 

! 

3591.3877258 

3.5 


1.1587874 

J 

4386.5308228 

3.6 


1.2601225 

! 

5357.7207642 

3.7 


1.3560715 

! 

6543-9348145 

3.8 


1.4477527 

! 

7992.7799683 

3.9 


1.5358592 

! 

9762.4033203 

4.0 


1.6206077 

! 

11923.8262939 
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